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ABSTRACT 

We perform linear analysis of thermal instability in a contracting large cloud filled with warm H I gas and 
investigate the effect of metallicity and radiation flux. When the cloud reaches critical density nt, the cloud 
fragments into cool, dense condensations because of thermal instability. For a lower metallicity gas cloud, the 
value of nf is high. Collision between condensations will produce self-gravitating clumps and stars thereafter. 
From the result of calculation, we suggest that high star formation efficiency and bound cluster formation are 
realized in low-metallicity and/or strong-radiation environments. 

Subject headings: galaxies: star clusters — globular clusters: general — instabilities 
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1. INTRODUCTION 

The Galactic globular clusters are suggested to be fossil 
record of the Galaxy formation era (e.g., Fall & Rees 1985). 
So any knowledge about their formation mechanism will be 
helpful to understand the protogalactic environment. Young 
counterparts of globular clusters are absent in the Milky Way, 
whereas in galaxies such as the LMC, SMC and M33, young 
globulars exist (e.g., Kumai, Basu, & Fujimoto 1993; Efremov 
1995). Such galaxies are in many case ongoing mergers and 
merger remnants or starburst galaxies, but young globulars also 
exist in normal spiral, dwarf and irregular galaxies (Schweizer 
1999). The conditions that enable globular cluster formation in 
these galaxies may be also realized at the epoch of the Galactic 
globular cluster formation. Since globular clusters are gravi- 
tationally bound, we attempt to clarify the condition of bound 
cluster formation. 

At the formation epoch of the globular clusters, high star for- 
mation efficiency (SFE) must be achieved somehow, since high 
SFE is required in order to form a bound cluster. The thresh- 
old SFE value for bound cluster formation is 0.5 in the case 
of rapid gas removal (e.g., Hills 1980), which is significantly 
greater than average SFE (~ 0.01) in the giant molecular clouds 
(GMCs). Thus, to understand the conditions of bound clus- 
ter formation, we should know in what environment the SFE 
can be high. In regard of this problem, Elmegreen & Efre- 
mov (1997) claimed that high-pressure environment is required 
for bound cluster formation because disruption of the cluster- 
forming cloud by massive stars is difficult in such an environ- 
ment. They also argued that bound cluster formation is efficient 
in a low-metallicity environment, as it is difficult to disrupt the 
cloud by radiation pressure of massive stars. Here, considering 
the formation process of dense clouds, we point out that high 
SFE and efficient bound cluster formation can be achieved in 
low-metallicity and/or strong-radiation environments because 
of cloud fragmentation via thermal instability. 

Molecular clouds, where stellar clusters form, usually be- 
long to some larger structure such as superclouds (Elmegreen 
& Elmegreen 1983), the largest cloud complexes in galaxies. 
Due to the similarity between the masses of superclouds and 
the critically unstable Jeans mass, large-scale gravitational in- 
stability is considered to produce these superclouds (Elmegreen 



1987). Warm H I gas occupies large fraction of these super- 
clouds, since roughly half of the H I in the Milky Way disk is 
in the form of warm gas (Kulkarni & Heiles 1988). When the 
gas disk becomes gravitationally unstable and fragments into 
superclouds, this warm gas should be converted into cold H I 
gas since the gas becomes thermally unstable when the density 
is somewhat higher. From this cold H I gas, GMCs and, subse- 
quently, stellar clusters are expected to form. Similarly, GMC 
formation in a supercloud via thermal instability was also con- 
sidered by Kolesnik (1991). However, his method was based on 
density structure of static isothermal cloud and our treatment is 
different. 

For modeling the above-mentioned top-down scenario of 
cluster formation from a supercloud (see also Efremov 1995), 
we consider a large, warm gas cloud (T ~ 10 4 K) contracting 
because of its self-gravity. As the cloud density increases, it 
becomes thermally unstable and breaks up into cool, dense con- 
densations. We investigate this process based on linear analysis 
of thermal instability in the cloud. For the evolution of conden- 
sations, we adopt the scenario that these condensations collide 
each other and self-gravitating clumps are formed in the cloud, 
and then stellar clusters form. To judge whether the cluster is 
bound or not, we should know the mean density of the cloud at 
the point of breakup, since we can relate the mean density of 
the cloud to its velocity dispersion, which determines SFE and 
hence whether the cluster is bound or unbound. So we calculate 
the thermal evolution of a contracting gas cloud and estimate 
the density when the cloud becomes thermally unstable. The 
effects of metallicity and radiation flux are also examined. 

In the following sections, we present the model of the cloud 
and the instability analysis method. Based on the results of the 
calculations, we discuss the conditions necessary for high SFE 
and bound cluster formation. 



2. THERMAL EVOLUTION OF A CONTRACTING GAS CLOUD 

We calculate the evolution of density, temperature, and ion- 
ization degree of a contracting gas cloud. For simplicity, we 
consider a spherical and homogeneous gas cloud. As an ini- 
tial condition, we take the static cloud filled with warm H I gas 
(T ~ 10 4 K) in thermal and ionization equilibrium. The con- 
traction is simplified as a free-fall collapse. Then the thermal 
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evolution of the cloud is given by the following basic equations: 
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Here p and n are the mass and number density of hydrogen of 
the cloud, respectively, and «, is the initial value of n. P, T, e, 
and x e are the pressure, temperature, internal energy, and ion- 
ization degree of the cloud, respectively. Function C represents 
the net cooling rate. As heating sources, we consider inter- 
stellar UV radiation, soft X-ray, and cosmic ray. As coolants, 
we consider hydrogen Lya cooling, line cooling by O I, C II, 
Si II, Fe II, dust cooling, and so on. Ionization states of these 
metal species are fixed according to their ionization potentials. 
In this work we approximate the cloud as optically thin. Func- 
tion X represents the net ionization rate of hydrogen atoms. We 
consider soft X-ray and cosmic ray as ionizing sources. To 
determine the X-ray heating and ionization rate, the hydrogen 
column density of the cloud, Nh, must be specified, we take 
A^h = 10 19 cm" 3 throughout the cloud evolution, though A^h in- 
creases as the cloud contracts. This simplification is not serious, 
since we performed the same calculation including the change 
of A^h and got almost the same result. For these thermal pro- 
cesses, we refer to Wolfire et al. (1995) and references there 
in. 

To verify the dependence on metallicity and radiation flux, 
we introduce parameters Z and Xq. In our assumptions, the 
number density of dust grains and the gas phase abundances of 
O I, C n, etc. change in proportion to Z. The UV radiation flux, 
X-ray flux and cosmic ray flux change proportional to Xq. Here 
Z = 1 and Xo = 1 correspond to solar values. 

3. FORMATION OF CONDENSATIONS DUE TO THERMAL 
INSTABILITY 

We estimate the mass and the formation timescale of con- 
densations using linear approximation. The growth of density 
perturbation in the parent cloud is governed by the equation 
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where k and a are the wave number and the growth rate of the 
perturbation, respectively. This equation is the same as in Field 
(1965), except in the last term, where it represents nonequilib- 
rium effect. Thus, the net cooling rate Co = C(To,po,x e o) is 
not necessary zero, where 7b, po, and x e Q are the background 
temperature, density, and ionization degree, respectively. The 
definitions of k„, k T , and k K are the same as in Field (1965): 
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where p, 7, R, and c s are the mean molecular weight, ratio of 
specific heats, gas constant, and sound speed, respectively. We 



take the value of p as 1.4 and 7 as 5/3. The sound speed is ex- 

pressed as c s = (7P0/P0) . where Po is the background pres- 
sure. The symbols C p , Cj, and Kq are represented as follows: 
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where K(T) is thermal conductivity, and we adopt K(T) = 
2.0 x 10 3 r 1//2 , the thermal conductivity for neutral hydrogen. 

Among three roots of equation (4), two roots correspond to 
sound wave modes and one root corresponds to isobaric con- 
densation mode (Field 1965). When this equation has unstable 
isobaric condensation mode, condensations form because of the 
growth of unstable perturbations (e.g., Goldsmith 1970). This 
mode has a maximum growth rate of a m at wave number k m 
(Fig. 1). Among various density perturbations with different 
intensities and wave numbers in the cloud, the perturbations 
with wave number k m grow fastest. So we estimate the mass of 
the condensation, M, as 

(7) 
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We also estimate the formation timescale as 

t g =—. (8) 

Because equation (4) is derived under the assumption that 
contraction of the cloud is negligible, it must keep in mind 
that equation (4) is applicable only when the growth timescale 
of perturbation is shorter than the contraction timescale of the 
cloud. Also note that our use of equation (4) is justified when 
A m is much smaller than the parent cloud, because equation (4) 
is for an infinite homogeneous medium. We can neglect this 
point, since the parent cloud (> a few kpc ) is much larger than 
\ m ( ~ a few pc). 

4. RESULTS 

We perform the calculation in the parameter ranges Z = 1- 
10" 2 , X = 10-0.1. For the purpose of illustration, we show 
the results of the cases (Z = l,X = 1) and (Z = 10" 2 ,X = 1) 
in Figures 2 and 3. In the figures, evolutionary paths in n— T 
plane of both the free-fall collapsing gas cloud and the quasi- 
statically contracting gas cloud keeping thermal equilibrium are 
presented. Equations (l)-(3) are used for the free-fall collapsing 
paths and the equations, C(T,p,x e ) = and X(T,p,x e ) = 0, are 
used for thermal equilibrium paths. The condensation mass M 
and tg/tdyn, the ratio of formation timescale of condensation to 
dynamical timescale of the parent cloud, are also shown in the 
same figures. Here the dynamical timescale of the parent cloud, 
fdyn, is defined as pj (dp/dt), where p is the density of the par- 
ent cloud. The contours of M and t g /u yn in Figures 2 and 3 are 
calculated using equation (4), with the assumption of ioniza- 
tion equilibrium (i.e., only n and T are used for calculation, and 
x e is a function of n and T). From the figures we can see the 
tendency that the smaller metallicity, the more the evolution- 
ary path of the free-fall collapsing cloud deviates from the path 
of the quasi-statically contracting cloud. For the case Z = 1, 
the two evolutionary paths deviate in temperature ~ 1000K at 
density ~ 1 cm" 3 . In contrast, for the case Z = 10~ 2 , the two 
paths deviate ~ 8000K at density ^100 cm" 3 . For the paths of 
the free-fall collapsing gas cloud, dependence on metallicity is 
stronger than equilibrium paths. Such behavior is explained by 
the fact that a nonequilibrium state is realized for free-fall cases 
because the cooling timescale can be long relative to dynami- 
cal time and that degree of nonequilibrium is stronger for lower 
metallicity since cooling time is in inverse proportion to Z. 
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FIG. 1 . — Dispersion relation of thermal instability (schematic example). Among three modes, only condensation mode is presented. Since the region a > 
exists, this mode is unstable and condensations are formed from density perturbations with positive a. Note that there exists a maximum growth rate a m at wave 
number k m . 



As the cloud density exceeds some critical value, the gas 
cloud becomes thermally unstable. At this point, the growth 
timescale of perturbation, t g , is longer than the dynamical 
timescale of the cloud, fdyn- So the perturbation cannot grow to 
the condensation, because the parent cloud evolves faster than 
the perturbation. As the cloud evolves, t g becomes smaller than 
fdyn- Thus, condensation will form at the moment t g = fd yn - For 
the case Z = 1, condensations form at the density ~ 1 cm" 3 (Fig. 
2a), whereas for the case Z = 10" 2 , condensations form at the 
density <~ 300 cm" 3 (Fig. 3a). Condensation mass is a few 
1OM at the density <~ 1 cm" 3 for Z = 1 (Fig. 2b) and a few M Q 
at the density ~ 300 cm" 3 for Z = 10" 2 (Fig 3b). 

For the case when Z and Xq are simultaneously changed, the 
results are shown in Figure 4. The calculational procedure for 
Figure 4 differs from Figures 2 and 3, in a point that t g and M 
are calculated using equations (l)-(8) to draw Figure 4, whereas 
the assumption of ionization equilibrium is adopted to draw the 
contours of Figures 2 and 3. In Figure 4a, the cloud density «/ 
at the time of condensation formation is shown. For Xq = 1, the 
density rif increases from ~ lcm" 3 at Z = 1 to <~ 300cm~ 3 at 
Z = 10~ 2 . For Z = 1, rif increases from <~ 0.1cm" 3 at Xq = 0.1 
to ~ 10cm" 3 at Xq = 10. In Figure 4b, the condensation mass is 
shown. ForXo = 1, the mass M decreases from <~ 30M Q at Z = 1 
to ~ IM Q at Z = 10" 2 . For Z = 1, M decreases from ~ 300M Q 
at Xq = 0.1 to ~ IM Q at Xq = 10. Thus, from Figure 4, we can 
see that for the smaller metallicity, the cloud density rif is the 
larger, and the condensation mass M is the smaller. We can also 
see that for the stronger radiation, the density rif is the larger 
and the mass M is the smaller. 

Previously, Elmegreen & Parravano (1994) calculated the 
critical pressure for coexistence of cool and warm gas, chang- 
ing radiation flux and metallicity. This critical pressure roughly 
corresponds to n / in this work. Contrary to our results, metal- 
licity dependence of the critical pressure is small in their calcu- 
lation. They assumed static cloud for the background, whereas 
we assume free-falling cloud. Then, as noted above, for the 
smaller metallicity, the evolutionary path of free-fall contract- 
ing cloud deviates more from the path of the quasi-statically 
contracting cloud, because of the reduced cooling rate. On the 
other hand, for the static cloud, thermal equilibrium is achieved 
so that metallicity dependence is weak. This is the reason of the 
difference between our work and their work. 

We estimate the density of condensation as follows. After 



the criterion t g = fd yn is achieved, the overdensity region (size 
~ A m ) in the parent cloud grows in pressure equilibrium with 
surrounding gas, since crossing time is short compared with 
cooling time for the fastest growing perturbation. Then the 
overdensity evolves along an isobaric path starting from the 
point t g = fdyn (this path is not drawn in Figs. 2 and 3). The 
runaway growth of the overdensity continues until it reaches 
thermochemical equilibrium, where it is thermally stable (e.g., 
Parravano 1987). We search the intersection of the isobaric path 
(the path with nT = constant, starting from the point t g = ?d yn ) 
with the equilibrium path (Figs.2 and 3, dashed line) and iden- 
tify it as "condensation". From the intersection, we can know 
the density of condensation, n c , and the temperature of the con- 
densation, T c . We perform the calculation and find that n c is 
roughly 100 times larger than rif in most parameter ranges. In 
the above-described procedure, we assume that the intersection 
is thermally stable and self-gravity of the overdensity can be 
neglected. We confirm that all intersections in this calculation 
are thermally stable, so that the first assumption is justified. The 
second assumption is also justified by the later discussion. Us- 
ing the values of n c and T c , the Jeans mass can be derived, and it 
is much larger than the condensation mass for most parameter 
values. Thus, self-gravity of the condensation is negligible and 
the density inside the condensation is expected to be almost uni- 
form. After formation, condensations accrete surrounding gas 
and can potentially grow to be self-gravitating. To estimate the 
mass growth rate of condensation, we should calculate the rate 
of mass flow across the boundary between the condensation and 
surrounding gas. According to previous works (e.g., Penston & 
Brown 1970; Parravano 1987), the growth timescale of the con- 
densation mass is much larger than the free-fall timescale of the 
parent cloud. So we can neglect the mass growth of condensa- 
tion through accretion. 

5. APPLICATION TO BOUND CLUSTER FORMATION 

After formation of condensations, the parent cloud contin- 
ues free-fall contraction until collisions between condensations 
begin. Then the formation of massive self-gravitating clumps 
through collisional buildup of small condensations is expected. 
Star formation inside such massive clumps will follow. How- 
ever, evaluation of such evolution is difficult task. The onset of 
collisions depends on velocity distribution and space density of 
condensations, which we cannot know precisely. Subsequent 
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FIG. 2. — (a) Ratio of formation timescale of condensation, t g , to dynamical timescale of the cloud, ta yn , is shown in n—T plane as solid contours. The dynamical 
time tiyn is defined as po/ (dpo/dt), where po is the density of the cloud. Numbers near the contours represent the values of t g /ti yn . Metallicity and radiation 
parameters are Z = 1 and Xq = 1, respectively. The thermally unstable region is inside the most outer contour. We calculate the contours under the assumption of 
ionization equilibrium. Thermal evolution of the free-falling cloud is presented in the plane as a dotted line. For reference, we also show thermal evolution of the 
quasi-statically contracting cloud keeping thermal equilibrium as a dashed line, (b) Condensation mass M(Mq) is shown in n-T plane as solid contours for the 
same parameters. Numbers near the contours represent the values of M(Mq). Meanings of dashed and dotted lines are the same as in Fig. 2a. 
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FIG. 4. — (a) Cloud density n/(crrT 3 ) at the time of condensation formation is presented in Xq-Z plane as contours. Numbers near the contours represent n/fcm" 3 ). 
We assume that condensations form at the moment t g = ti yn , where t g is the formation timescale of condensation and «d y n is the contraction timescale of the cloud, 
(b) Condensation mass M(M Q ) at the time of its formation is presented in X -Z plane as contours. Numbers near the contours represent M(Mq). 



evolution of condensations also depends on velocity and mass 
distribution of condensations through the collisional process. 
This collisional process is not so simple (e.g., Hausman 1981) 
and makes the calculation of evolution complicated. In regard 
of these difficulties, we adopt a simplified picture of the evolu- 
tion of condensations. We neglect the random velocity of con- 
densations and assume homologous contraction of the parent 
cloud composed of condensations. Then the collisions begin 
when the mean parent cloud density increases to the density of 
condensation, n c , because the parent cloud is almost filled with 
condensations in such situation. As a result of collisions (we as- 
sume that the output of collision is mostly adhesion), massive 
self-gravitating clumps will form. At this point, the density of 
massive clumps is larger than the mean parent cloud density n c , 
because of self-gravity. So the parent cloud contract further. 
We assume collisionless collapse of the system composed of 
self-gravitating clumps, for simplicity. Then the virialization of 
the system and the beginning of star formation in the clumps is 
expected. 

It is reasonable to consider that the cluster-forming region is 
only a part of the parent cloud, not the entire cloud. We de- 
fine the mass of the region as M Q \ and take the reference value 
as 1O 5 M according to typical mass of globular clusters. Then 
virial velocity of the region is 

/ M cl \ 1/3 / n d \ [ / 6 , 

c ~ 2x {wk) (1^) kms < (9) 

where n z \ is the mean hydrogen number density of the cluster- 
forming region. Because (1) we assume collisionless collapse 
of the parent cloud composed of self-gravitating clumps, where 
the mean cloud density at the onset of collapse is the same as 
condensation density n c , and (2) condensation density n c is typ- 
ically 100 times higher than the parent cloud density at the time 
of condensation formation, n / (see the last paragraph in § 4), 
n c i and rif are related as 

« c i-800n/, (10) 



where the relations, n c \ ~ 8n e from the point (1) and n c ~ 100n/ 
from the point (2), are used. 

Whether the cluster is bound or not depends on SFE, the ra- 
tio of stellar mass to the mass of cluster-forming region. And 
SFE is determined by the interaction between stellar activity 
and the cluster-forming region. For example, stellar wind, radi- 
ation pressure, ionization pressure of the H II region or super- 
nova explosion may blow off remaining gas and lead to lower 
SFE. In this paper we consider two cases, gas removal by the 
expanding H II region and supernova explosion. 

First we discuss gas removal by expanding H II region. The 
pressure of H II region corresponds to velocity dispersion cuu ~ 
10km s" 1 . If virial velocity of the region, c, is smaller than chh. 
the H II region pushes out the surrounding gas and blows it off, 
and then the cluster will be unbound. In opposite case, the H II 
region cannot blow off the surrounding gas and high SFE can 
be achieved; then the cluster will be bound. So the criterion of 
bound cluster formation is considered as c > c H n. This corre- 
sponds to the condition 

B/ > 10 x f-g-V' ( CH " ^Y^V^m- 3 , (ID 

f vsooy Viokms-v vio 5 M y ' v ; 

where C represents density enhancement, i.e., the ratio of n c i to 
rif. Note that derivation of equation (11) relies on some rough 
assumptions and depends on uncertain factors such as C, Chii, 
and Ma- Moreover, dependence on cuu is very strong. So we 
should regard equation (1 1) as qualitative criterion. 

Next we discuss gas removal by supernova (SN) explosion. 
The proto-globular cluster cloud is disrupted by several to tens 
of concurrent SN explosions (Dopita & Smith 1986; Morgan & 
Lake 1989). A high star formation rate at the formation epoch 
of globular cluster will cause such multiple SN explosions so 
that relic gas will be blown off. When the star formation time 
scale ?sf is shorter than the lifetime of massive star, r, a large 
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fraction of gas in the cluster-forming region can be converted 
into stars before SN explosions; then high SFE is achieved and 
a bound cluster can be formed. In contrast, when tsF is longer 
than t, SN explosions will remove surrounding gas and stop 
later star formation. In this case, low SFE and unbound cluster 
formation is expected (e.g., Yoshii and Arimoto 1987). The star 
formation timescale fsF is expected as 

ts¥=At ff , (12) 

where A represents a proportional constant of order 1 and 

1 /2 

t/f = (3ir/32Gp) is free-fall timescale of the cluster-forming 
region. We adopt the lifetime of massive star, t, as 

r~5xl0 6 yr- (13) 

The star formation timescale fsF, evaluated at density n c \, is 

t SF =At ff ~ 1.3 x 10 7 A (_g-3)" 1/2 yr. (14) 
Thus, the condition for high SFE (r > ?sf) in term of rif is 

^ >10 K^kr) 2 (l) 2 (^) lcm " 3 - (15) 
This criterion is qualitative, since it depends some uncertain 
factors, the same as equation (11). 

For both case, rif > 10 cm" 3 is required for high SFE and 
bound cluster formation. From Figure 4a we can see that the 
density n f increases as metallicity decreases and also as radia- 
tion increases. The condition rif > 10 cm" 3 corresponds to the 
parameter range Z < 0.05 for Xq = 1 and Xq > 8 for Z = 1 . This 
means high SFE and efficient bound cluster formation are re- 
alized in a low-metallicity and/or strong-radiation environment 
in our model. From Figure 4a, we can also see that when Z ~ 1, 
n/ is more sensitive to Xq than Z and when Z < 0.1, rif is more 
sensitive to Z than Xq. Thus SFE mainly depends on Xq around 
Z ~ 1 and on Z below Z ~ 0. 1 . 

The qualitative behavior of the above result will be inde- 
pendent of the details about evolution of condensations and 
feedback from star formation, since the essential ingredients 
are only two points, (1) low-metallicity and/or strong -radiation 
fields lead to high breakup density of the parent cloud, and (2) 
high-density environments result in high SFE. However, our 
conclusion may be changed if we consider nonspherical evo- 
lution of the cloud. The dynamical timescale of an initially 
nonspherical cloud or a cloud with rotation is similar to that of 
a spherical cloud at first, e.g., fdyn °c n" 1 / 2 , where n is the den- 
sity of the cloud. But later disklike collapse will follow and 
the dynamical timescale becomes fd yn oc n" 1 (Susa, Uehara, & 
Nishi 1996). Then the ratio of the cooling time t coo \ oc n~ l Z~ l to 
the dynamical time is 

^COOl -1/2 ry-l /, 

oc n ' 27 , (16) 

fdyn 

for initial spherical-like phase and 

^ocZ" 1 , (17) 

^dyn 

for later disklike phase. When the cloud becomes thermally 
unstable, we can estimate whether or not condensations form, 
using the value of / C ooiAdyn- If this ratio is less than unity, con- 
densation will form, and vice versa. In case the cloud becomes 
unstable when it is still in a spherical-like phase, condensations 
will form regardless of the metallicity Z, as we calculated in 



this work. This can be explained since the density n, when 
the cloud becomes unstable, increases as Z decreases (see Figs. 
2 and 3), so that / CO oi/%n does not change so much. But, in 
case the cloud becomes unstable when it is in a disklike phase, 
condensations may not form for the low-metallicity case since 
fcooiAdyn increases as Z decreases. In this case our discussion 
in this section is not applicable. To summarize, in the presence 
of rotation or deviation from spherical symmetry, if the cloud 
becomes thermally unstable in the early phase of collapse, con- 
densations will form and our conclusion in this section is ade- 
quate. In contrast, if the cloud becomes thermally unstable in 
the late phase of collapse, condensations may not form through- 
out the collapse and another discussion may be required. 

Finally, we present some applications. A low-metallicity en- 
vironment is realized in dwarf galaxies and in the early stage of 
our Galaxy, for example. And strong-radiation fields are real- 
ized in starburst galaxies and probably in the early stage of our 
Galaxy, where young bulge stars are expected to provide strong 
radiation. Thus, according to the result, young globular clus- 
ters in dwarf galaxies and starburst galaxies may be produced 
through the mechanism we show in this work. Also, the disk 
population of globular clusters in our Galaxy (Zinn 1985) may 
be produced through the same mechanism. 

6. SUMMARY 

In this paper we perform linear analysis of thermal insta- 
bility in a collapsing warm gas cloud and study the effect of 
metallicity and radiation flux. For the calculation, we adopt 
a one-zone approximation and the cloud contraction is simpli- 
fied as free-fall collapse. When the cloud density reaches crit- 
ical value rif, the cloud fragments into cool dense condensa- 
tions via thermal instability. From our calculation, the critical 
density n / increases as metallicity decreases, and also as radia- 
tion increases. Condensations collide with each other and self- 
gravitating clumps will be produced when the mean cloud den- 
sity becomes sufficiently high; then stars will form. Expansion 
of the H II region around the massive star and supernova ex- 
plosions will blow off surrounding gas and stop star formation 
process. When the mean density at the time of star formation 
is high, high virial velocity prevents expansion of the H II re- 
gion. Also, in such high-density environments, the star forma- 
tion timescale is shorter than the lifetime of a massive star. Then 
the gas in cluster-forming region will be converted into stars ef- 
ficiently, before the gas is dispersed by expanding H II region or 
supernova explosions. In our calculation, such high density is 
realized in the contracting low-metallicity gas. And if the for- 
mation of a contracting gas cloud is possible, a strong-radiation 
environment is another candidate. Thus, it is suggested that 
high star formation efficiency and bound cluster formation are 
likely achieved in low-metallicity and/or strong-radiation en- 
vironments. Such environments exist in dwarf galaxies, the 
early stage of our Galaxy and starburst galaxies. According to 
the result, globular clusters in these galaxies may be produced 
through the mechanism we show in this paper. 
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